function main()
% MAIN - 计算程序
%   
    % Watson 函数
    c=clock;
    b=c(4)*3600+c(5)*60+c(6);b=round(b*10);
    filename=num2str(b);
    filename=['result/',filename,'.result'];
    fid = fopen(filename, 'w');

    N = [2 10 20 31];
    %    N = [2 6 9 12];
    fprintf(fid, 'Watson:\n');
    fprintf(fid, '%6.0f',N);
    fprintf(fid, '\n');
    [F,info]=WatsonTest(N);
    
    fprintf(fid, '\t\tdamping\t\tmulti\t\tSR1\t\tBFPS\t\tDFP\n');F=F';
    fprintf(fid,  '\t\t%12.8e\t\t%12.8e\t\t%12.8e\t\t%12.8e\t\t%12.8e\n',F);
    fprintf(fid,'\n\n');
    for i = 1:length(N);
        out = info(:,:,i)';
        fprintf(fid,  '%6.0f %6.0f %6.0f %6.0f %6.0f\n',out);
        fprintf(fid,'\n\n');
    end
    
    % Discrete Boundary Value 函数,
    N = [4 20 40 60];
    fprintf(fid, 'Discrete Boundary Value 函数:\n');
    fprintf(fid, '%6.0f',N);
    fprintf(fid, '\n');
    [F,info]=DBVTest(N);
    fprintf(fid, '\t\tdamping\t\tmulti\t\tSR1\t\tBFPS\t\tDFP\n');F=F';
    fprintf(fid,  '\t\t%12.8e\t\t%12.8e\t\t%12.8e\t\t%12.8e\t\t%12.8e\n',F);
    fprintf(fid,'\n\n');
    for i = 1:length(N);
        out = info(:,:,i)';
        fprintf(fid,  '\t\t%6.0f\t\t%6.0f\t\t%6.0f\t\t%6.0f\t\t%6.0f\n',out);
        fprintf(fid,'\n\n');
    end
    
    
    %Biggs EXP6 函数
    
    N = [6 10 20 30 50];
    fprintf(fid, 'Biggs EXP6 函数\n');
    fprintf(fid, '%6.0f',N);
    fprintf(fid, '\n');
    [F,info]=Biggs6Test(N);
    fprintf(fid, '\t\tdamping\t\tmulti\t\tSR1\t\tBFPS\t\tDFP\n');F=F';
    fprintf(fid,  '\t\t%12.8e\t\t%12.8e\t\t%12.8e\t\t%12.8e\t\t%12.8e\n',F);
    fprintf(fid,'\n\n');
    for i = 1:length(N);
        out = info(:,:,i)';
        fprintf(fid,  '%6.0f %6.0f %6.0f %6.0f %6.0f\n',out);
        fprintf(fid,'\n\n');
    end
    fclose(fid);
end
